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Abstract 

A new algorithm for extraction of asymmetries from polarized lepton- 
nucleon scattering data is proposed. The algorithm is stable to set-up ac- 
ceptance and/or luminosity monitor acceptance variations. A statistical test 
for checking the data quality is proposed. 
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1 Introduction 
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The measurement of asymmetry in deep inelastic scattering of polarized lep- 
tons on polarized nucleons gives the main information for studying the internal 
spin structure of nucleons. This approach allows to avoid the problems of nor- 
malizations and, to some extend, acceptance corrections to obtain then polarized 
structure functions. 

Usually asymmetry is evaluated according to the formula: 

1 N+-N- 

A = (I) 

pt P b N+ + N- ' 1 ; 

where N + (N~) is the number of scattered leptons for the target spin parallel (anti- 
parallel) to the beam spin orientation, P l is a value of target polarization and P b 
is a value of beam polarization. This formula is valid for the case when luminosity 
and the set-up acceptance for the parallel and anti-parallel spin orientation is the 
same [Q]. Usually this condition takes place when two targets with the opposite 
spin orientation are used for measurements . 

A more general formula is used in || where measurements are carried out with 
only one target, sequentially for each target spin state: 

N-L+ - N+L~ 
N-L% + N+L P 

where L + (L~) is luminosity for each target spin state, and Lp(Lp) is luminosity 
weighted by the product of the beam and target polarization values. The for- 
mula (2) is valid if the set-up acceptance is the same and the luminosity monitor 
acceptance is the same for all states of the target spin. 

Another method based on exact formulae for differences of cross-sections for 
parallel and antiparallel spin orientations is developed in ||. According to this 
approach the number of events Nij collected for any kinematical bin and for a 
given pair (bin) of the beam (i) and target (j) polarizations, can be written as 

= a u L tj + crpLijP^ , (3) 

where Lij is the luminosity, P b , (Pj) is the measured beam(target) polarization, 
a u is the unpolarized part of the cross section and a v is the polarized part of the 
cross section. Acceptance, dead time and other corrections are assumed to be 
known and properly accounted for. The unpolarized part of cross section a u and 
the polarized part of cross section a p are extracted from the relationship (3) by 
using the likelihood method. Asymmetry is calculated as: 

A=^ . (4) 

In practical cases the set-up and luminosity monitor acceptances are not known 
and they vary with time. It can be the cause of the asymmetry bias if the methods 
[01, IHI ase used. 
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In this paper, we propose an algorithm for asymmetry extractions that is stable 
to acceptances variations. In framework of this algorithm a statistical test of 
the data quality check is also developed. The usage of the algorithm does not 
assume that acceptances are known or acceptances does not vary with time when 
measurements are performed. 

The paper consists of 4 sections. In Sect. 2 we describe an algorithm for 
asymmetry extraction. An example of using the algorithm is described in Sect. 3. 
Conclusions are given in Sect. 4. 

2 Algorithm for an asymmetry extraction 



Usually the data to determine asymmetry in lepton-nucleon interactions are 
collected during long periods of time , it can be done within several years. The 
beam and target polarization, beam intensity, luminosity, spectrometer accep- 
tance, luminosity monitor acceptance, and other experimental conditions vary in 
these periods of time. To minimize the systematic errors associated with these 
variations, the data are divided in samples with relatively stable experimental 
conditions and the target polarization is frequently reversed. 

The number of events iV collected for any kinematical bin for each sample can 
be written as: 

N = + (5) 

c c 

where C is the set-up acceptance, c is the luminosity monitor acceptance, I is 
the measured luminosity and P*(P b ) target (beam) polarizations. The measured 
luminosity can be presented as: 

1 = Lc , (6) 

here L is an unknown true value of the luminosity. For simplicity we have assumed 
that the dead time corrections are known and properly accounted for. We have 
also assumed that the target is mononuclear. 

Let us define a supersample as a subset of samples following one by one in 
time. We also assume that for this subset the time dependence of the variable: 

s = , (7) 

c 

can be approximated by a line 

s(t) = a + bt (8) 

where, a and b, are free parameters. 

Assume that the set of samples can be divided in n supersamples and each i-th 
supersample, i = 1, n contains mi samples. The number of events iVY,- collected 
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for any kinematical bin and j'-th, j=l,...,rrii, sample of the z-th supersample, can 
be written as: 

N id = {at + biti^kj + A{ ai + h;i; i )l; J Pl J l% . (9) 

Our aim is to extract asymmetry from (9) taking into account statistical errors 
of Nij, luminosity, and polarizations. The solution of this problem is found in 
the framework of the maximum likelihood method. The task of the maximum 
lakelihood procedure is to find the asymmetry, true values of luminosities <kj>, 
target polarizations <P i *>, beam polarizations <P;j> and parameters cij, hi. 

If fluctuations of polarizations and luminosity are assumed to be Gaussian 2 and 
for to be Poisson, the logarithm of the likelihood functional L, is as follows: 



- InL = X! <kj> ( a i + htij) + A < l ij><Ptj><Ptj> ( a i + hUj) - 
i=i j=i 

-.V, ; /// (</,,:.' (a; + bitij) + A <Z y ><fJ.><fg.> (a* + + 

where 5Pj*-, are statistical errors of Zy, P^-, P^ . 

Minimization of Eq. (10) gives asymmetry A, and all other parameters <P/->, 
<Pij>, <kj>, cii,bi. The complete matrix of statistical errors is calculated as 
the inverse matrix of the second derivative matrix of the functional (10) at its 
minimum Q. 

Residuals examination can be applied to the quality check of the data used for 
asymmetry extraction ||. They are defined as: 

Res Ny = Nii J"? , (11) 

where the fitted value of the number of events 

iV-f =<kj> {aj + kUj) - A <kj><Pij><PPj> (en + htij) (12) 



and 



Pes P 



pt._ <p*.> 



spt. 



2 There is a tendency for errors that occur in many real situations to be normally distributed 
due to the Central Limit Theorem. If an error is a sum of errors from several sources, then no 
matter what the probability distribution of the separate errors may be, their sum will have a 
distribution that will tend more and more to the normal distribution as the number of components 
increases, by the Central Limit Theorem. Thus, the assumption of normality is not unreasonable 
in most cases. In any case we can later check the assumption by examining residuals. 
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The residuals plotted in the time order fluctuate around zero with the same vari- 
ation and have no time trend if the quality of the data is good and choice of 
supersamples is reasonable. 

3 Example of the algorithm usage 

In order to demonstrate the usage of the algorithm, pseudodata were simulated. 
We take asymmetry A equal to 0.1, a u — 1. , number of samples n equal to 80, 
luminosity Zj, i — 1, ...,80, polarizations Pf, P^, acceptances Ci and q as they are 
shown in Fig. 1, Using formulae (5) and (6), iVj and U are calculated. 

The measured number of events iV, is simulated as Poisson distributed random 
variable with a mean value equal to iVj. The measured values of luminosity li and 
polarization Pf, PP were obtained by adding fluctuations to the luminosity li and 
polarization P/, P\ values. Fluctuations were simulated as Gauss distributed ran- 
dom variables with the mean value equal to zero and a equal to 2% for luminosity, 
0.5% for target polarization and 5% for beam polarization. These pseudodata are 
shown in Fig. 2. 

The algorithm was used for processing the pseudodata divided into a different 
number of supersamples with the same number of samples. The result of the 
algorithm application is shown in Table 1. The last row of Table 1 presents the 
asymmetry extracted with the algorithm from M. 

In order to calculate the average value of the asymmetry and bias, the pseu- 
dodata described above were simulated 25000 times. The average value for the 
asymmetry <A> and bias defined as the difference between the asymmetry aver- 
age value and its true value 0.1 are also shown in Table 1. As follows from Table 
1., the algorithm proposed in this article gives an essentially lower bias than the 
algorithm in 

Residuals were analysed to choose the optimal number of supersamples. This 
analysis has shown that the number of supersamples can be chosen equal to 8 , 
residuals for this case are shown in Fig. 3. iVj and its fitted value N[ u are also 
shown in Fig. 3. 
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TABLE 1. Results of the asymmetry extraction for different numbers of super- 
samples. The x 2 are calculated for the number of events iVj and their fitted value 
iVf *. The last row presents the results of the algorithm |2j application. 



Number of 
supersamples 


A 


x 1 


<A> 


Bias 


1 


0.0799 ±0.1426 


89.4987 


0.1015 ±0.0009 


1.5% 


2 


0.0810 ±0.1426 


88.8028 


0.1016 ±0.0009 


1.6% 


4 


0.0679 ±0.1430 


81.8742 


0.0995 ±0.0009 


-0.5% 


8 


0.0584 ±0.1442 


63.0929 


0.1008 ±0.0009 


0.8% 


16 


0.0676 ±0.1451 


38.4059 


0.1007 ±0.0009 


0.7% 


Algorithm 


0.1112 ±0.1456 




0.1160 ±0.0009 


16.0% 



4 Conclusions 

The algorithm for extraction asymmetry from the polarized lepton-nucleon 
scattering data stable to the set-up and luminosity monitor acceptances is pro- 
posed. In framework of the developed algorithm a statistical test for checking the 
data quality is proposed. The algorithm has no restrictions related with small 
statistics of experimental events and can be applied for inclusive as well as for 
semi-inclusive physics. Principles of constructing the algorithm can be used in any 
other field of the experimental particle and nuclear physics. 
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Fig.l. The data for simulation of pseudodata: (a) luminosity, (b) target po- 
larization, (c) beam polarization, (d) set-up acceptance, (i) luminosity monitor 
acceptance. 
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Fig. 2. The pseudodata: (a) measured number of events, (b) measured luminosity, 
(c) measured target polarization, (d) measured beam polarization. 
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Fig. 3. Residuals (a) for number of events, (b) for luminosity, (c) for target 
polarizations, (d) for beam polarizations. 
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Fig. 4. The number of events for samples iVj and their fitted values N{ lt (close 
cycles). 



9 



